function p=P(y)
p=sum(abs(y).^2)/length(y);
end